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Solid-state nanopores are single molecule sensors that measure changes in ionic current as charged 
polymers such as DNA pass through. Here, we present comprehensive experiments on the length, 
voltage and salt dependence of the frequency of double-stranded DNA translocations through conical 
quartz nanopores with mean opening diameter 15 nm. We observe an entropic barrier limited, 
length dependent translocation frequency at 4M LiCl salt concentration and a drift-dominated, 
length independent translocation frequency at 1M KC1 salt concentration. These observations are 
described by a unifying convection-diffusion equation which includes the contribution of an entropic 
barrier for polymer entry. 


I. INTRODUCTION 

Polymer translocation through a narrow pore is a ubiqui¬ 
tous process in organisms. Classic examples include the 
passage of mRNA across the nuclear pore complex and 
DNA ejection by a bacteriophage. In vitro experiments 
on polymer translocation are possible using nanopores 
fabricated in solid-state materials or with reconstituted 
membrane proteins. The basic method for sensing the 
translocation of a polymer, such as DNA, relies on apply¬ 
ing a voltage across the nanopore and measuring changes 
in ionic current as molecules pass through. This simple 
premise underlies research for developing nanopores as 
biosensors for instance in next-generation sequencing ap¬ 
plications pP. Yet experimental data on polymer translo¬ 
cation phenomena present many puzzles that remain to 
be fully understood [2]. 

Conceptually, the electrophoretically driven transloca¬ 
tion of DNA through a nanopore consists of essentially 
three steps; (i) the DNA moves towards the pore entrance 
by a combination of diffusion and drift due to the elec¬ 
tric field outside the pore, (ii) the DNA is captured at the 
mouth of the pore, and (iii) the DNA threads through the 
nanopore causing a detectable ionic current change. All 
three steps taken together determine the translocation 
frequency defined here as the number of DNA strands 
passing through the nanopore per second per unit con¬ 
centration. Accurate measurements of the translocation 
frequency enable characterisation of the transport pro¬ 
cess and are important for applications in nucleic acid 
sensing. 

For the a-hemolysin nanopore, the translocation fre¬ 
quency of short, single-stranded (ss)DNA shows an expo¬ 
nential dependence on voltage under typical experimen¬ 
tal conditions indicating a barrier-limited process EHg. 
At high voltages and polymer concentrations, a second 
exponential is observed which is attributed to the effect 
of polymer-polymer interactions at the pore mouth [4]. 
Small diameter (sub 5-nm) solid-state nanopores show an 
increasing translocation frequency with double-stranded 
(ds)DNA length from 0.4 kbp to ~ 8 kbp followed by 
an indication of a length independent translocation fre¬ 
quency for longer lengths (6j. Electro-osmotic flow ef¬ 


fects on the dsDNA trapped at the membrane surface 
were suggested as a reason for the length dependence 
of translocation frequency 1313. It is known that the 
small diameter of these nanopores relative to the dsDNA 
cross-section results in strong interactions between the 
dsDNA and pore surface IS M- Larger nanopores, with 
diameters several times the cross section of dsDNA, show 
significantly less spread in transit times due to the lack 
of strong surface interactions J5]. This enables a sim¬ 
pler physical picture with the translocation speed accu¬ 
rately modelled by a combination of hydrodynamics and 
continuum electrostatics jlOj [IT]. While many reports 
in the literature have characterised the threading speed 
and conformations of dsDNA in these wide nanopores 
mm there are few measurements on the translocation 
frequency which can provide insight on the underlying 
transport mechanisms. 

In this paper, towards developing a complete picture of 
dsDNA translocation in solid-state nanopores, we study 
the voltage, length, salt and concentration dependence 
of dsDNA translocation frequency for nanopores with di¬ 
ameters 15±3 nm. We show the simultaneous detection 
of ten different DNA lengths in the range 0.5 kbp to 
10 kbp which allows us to accurately characterize the 
translocation frequency. In 4M LiCl electrolyte, there is 
an increase in translocation frequency with DNA length 
whereas experiments using 1M KC1 show a length inde¬ 
pendent translocation frequency. These trends are ac¬ 
curately captured by a ID convection-diffusion equation 
which incorporates an entropic barrier for polymer entry 
into the nanopore and scaling laws for polymer diffusion 
and mobility. 

II. METHODS 

Quartz capillaries were fabricated into nanopores using 
a previously published protocol with diameters 15 ± 3 
nm (meanis.d.) [151 . The nanopore was filled with a 
solution of 10 mM Tris-HCl pH 8, 1 mM EDTA and elec¬ 
trolyte (either 4M LiCl or 1M KC1). DNA was added to 
the reservoir outside the nanopore tip which was set as 
the electrical ground. Table I gives the final reservoir con¬ 
centrations of each DNA strand used for the experiments 
in 4M LiCl in Figures 2 and 4. The sample was made 
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DNA 

length (kbp) 

0.5 

1 

1.5 

2 

3 

4 

5 

6 

8 

10 

Cone (nM) 

3.2 

1.6 

0.91 

0.91 

1.6 

0.31 

0.32 

0.32 

0.20 

0.16 


TABLE I. Final reservoir concentration of DNA lengths in 
ladder used for Figures 2 and 4. 

by diluting a mix of DNA lengths (1 kb DNA ladder) 
purchased from New England Biolabs. For experiments 
using 1M KC1, three DNA lengths of 5 kbp, 10 kbp and 
20 kbp (NoLimits chromatography-purified DNA, Ther- 
moFisher Scientific) were mixed in an equimolar ratio so 
that the final reservoir concentration was 0.12 nM of each 
length. 

Fig. la shows a typical scanning electron microscope 
image of a conical quartz nanopore and in Fig. lb the 
electric field profile at +600 mV applied potential is 
shown from finite-element calculations. The electric field 
was calculated using a finite-element solver (COMSOL 
4.4) and the full Poisson-Nernst-Planck equations for ion 
transport through the pore in a 2D axisymmetric geom¬ 
etry. The electric field is mainly confined to the last few 
hundred nanometres of the tip and decays rapidly away 
from the pore entrance. Ionic currents were measured us¬ 
ing an Axopatch 200B amplifier, filtered at 50 kHz using 
an external 8-pole Bessel filter and subsequently sampled 
at 250 kHz. All data was analysed with a custom-written 
program in Labview 2013 (National Instruments). All ex¬ 
periments were performed at 20° C. Translocations were 
detected using a threshold analysis with minimum cur¬ 
rent change set to 50 pA, minimum translocation dura¬ 
tion 50 /rs and minimum event charge deficit of 3 fC. The 
translocations from all DNA lengths at all voltages fall 
outside these limits and therefore we do not miss translo¬ 
cations (Figure 15). 

III. RESULTS 

Initially we characterised the translocation of individual 
dsDNA lengths in 4M LiCl electrolyte. We observe that 
dsDNA can translocate in folded configurations in agree¬ 
ment with previous studies using similar sized nanopores 
mm- This is shown by the quantised levels of current 
blockade as for example in Fig. Id. For 10 kbp DNA we 
observe folded configurations with up to three double¬ 
strands of the same molecule in the pore (Figure 7). The 
translocation time is determined by both the DNA length 
and the folding configuration of the DNA and does not 
uniquely identify a single DNA length. However, the 
integral of the current with respect to time or event 
charge deficit (ECD), can be used to distinguish each 
DNA length. Fig. le shows ECD distributions of indi¬ 
vidual dsDNA lengths of 1 kbp, 5 kbp and 20 kbp which 
are single peaks each well fitted by a Gaussian function 
(Fig. le). Outliers are observed at low ECD values for 
both 5 kbp and 20 kbp which we attribute to a small 
amount of fragmentation from DNA shearing which can 
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FIG. 1. (a) Scanning electron microscope image of a typical 

quartz nanopore, (b) Finite-element analysis of the electric 
field profile at a voltage of +600 mV using the geometric 
model given in [I+j. Lines represent tangents to the electric 
field, (c) Typical current trace showing downward spikes due 
to DNA translocation interrupting the baseline. The event 
in the dashed area is shown at greater time resolution in (d) 
which also shows the event charge deficit (ECD) in blue - 
the integrated area of the event relative to the baseline, (e) 
Histograms of ECD for 1 kbp, 5 kbp and 20 kbp DNA with 
N=1119, N=764 and N=456 translocations respectively. Each 
length was measured separately at +600 mV with a DNA 
concentration of 1.6 nM. 

occur during pipetting of long DNA molecules m and 
has been observed previously IT2: . 

Our observation of a single peak in ECD for each 
DNA length strongly suggests that we only measure 
a current signal above the background noise when the 
molecule fully translocates through the nanopore. It 
is also possible that there are events where the DNA 
is captured at the pore mouth but does not overcome 
the energy barrier to translocation and therefore returns 
to the sample reservoir under thermal motion. These 
collisions can create a second population and have been 
observed for dsDNA in sub-5 nm solid-state nanopores 
mu, in a pressure-voltage trap system for 10 nm diam¬ 
eter nanopores [18| and for ssDNA in a-hemolysin nnj- 
We attribute the lack of a second population from colli¬ 
sions as being due to the conical shape of the nanopores 
meaning that the effective sensing length is on the order 
of a few hundred nanometres. Therefore collisions do not 
cause a significant enough current change to be measured 
above the noise. 

For accurate, simultaneous determination of transloca¬ 
tion frequency of a range of lengths, we used a mixture of 
10 DNA lengths in the range 0.5 kbp to 10 kbp (Fig. 2a). 
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FIG. 2. (a) Agarose gel (1%) electrophoretogram of the DNA 
sample containing 10 lengths in the range 0.5 kbp to 10 kbp. 
(c) ECD histogram for 5711 translocation events recorded at 
+600 mV. The histogram is logarithmically binned and the 
peaks corresponding to the shortest and longest DNA lengths 
are labelled, (c) Scatter plot of mean event current against 
event duration for same events as in (b). The dashed grey 
lines are fits according to AlMean = — ac ^ g ° l; and the solid 
line shows the 3 fC threshold. 


The ECD histogram for 5711 translocations with 4M LiCl 
electrolyte is shown in Fig. 2b and Fig. 2c shows the cor¬ 
responding scatter plot of mean current versus duration 
for each translocation. The low noise of quartz nanopores 
combined with the slow DNA velocity in 4M LiCl means 
that we have sufficient resolution to identify the 10 DNA 
lengths present with each DNA length in the scatter plot 
observed as a band. The dashed grey lines are given by 
AlMean = — ECI ^' !ak where AlMean is the mean event 
current, At is the event duration and ECDp ea k is the 
centre of the Gaussian function determined by fitting 
each ECD peak. The excellent fit of this equation to 
the data confirms that for all translocations of a particu¬ 
lar length, the ECD value is constant irrespective of the 
folding state of the DNA as it passed through. As de¬ 
scribed earlier, we assign each data point to a transloca¬ 
tion of DNA and therefore determined the total number 
of translocations for each length from the integrated area 
of Gaussian functions fitted to each peak in the ECD his¬ 
togram. The translocation frequency (shown in Fig. 4) 
was then calculated as N/ct where N is the number of 
translocations, c is the DNA concentration and r is the 
experiment time. 

IV. TRANSLOCATION FREQUENCY MODEL 

We use the following model to describe the DNA translo¬ 
cation frequency based on Muthukumar [20]. There are 
three basic steps for the translocation process. In the 


first step, the DNA molecule undergoes drift and diffu¬ 
sion as a whole molecule with its allowed conformational 
fluctuations towards the mouth of the nanopore where 
it gets caught. At the end of this capture stage, the 
effective diameter (twice the radius of gyration) of the 
DNA molecule is much bigger than the 15 nm mean di¬ 
ameter of the nanopore. Therefore, in the next step of 
successful translocation, the DNA chain can put only a 
few segments and not the whole chain. This step (second 
step) is entropically unfavorable due to the restriction 
of the number of conformations of a translocating DNA 
molecule. Hence the DNA molecule suffers from a free 
energy barrier. Also, since the pore diameter is larger 
than the thickness of DNA (2 nm), but not much bigger, 
several DNA segments can penetrate the mouth of the 
pore simultaneously. Depending on the particular man¬ 
ner in which these few segments are placed at the pore 
mouth, the values of the corresponding entropic barriers 
can vary mm- Since several confined conformations of 
DNA are allowed, unlike in single-file translocation, it is 
reasonable to assume an effective free energy barrier for 
the whole chain. 

It would be of course desirable to perform simulations 
to compute the particular values of the free energy barri¬ 
ers for various particular configurations of tails versus 
loops and multiple loops of DNA at the pore mouth. 
Since this is much beyond the scope of the present status 
of the subject, we approximate the net effect in an ensem¬ 
ble of translocation events as an effective barrier. Once 
this second step of nucleation is successful, the third step 
of threading away the molecule is down-hill in free en¬ 
ergy. In view of these considerations, we model a stochas¬ 
tic one-dimensional process of the polymer negotiating a 
free energy landscape where the first step is governed 
by diffusion and electrophoretic drift in bulk, the second 
step by an effective entropic barrier per chain, and the 
third step by the driving electric field. The net flux of 
molecules is determined by the entropic barrier, diffusion 
of the DNA and the electrophoretic force acting on the 
DNA. We can neglect electro-osmotic flows since mea¬ 
surements using optical tweezers show that the flow rate 
decreases strongly with salt concentration and is insignif¬ 
icant compared to the electrophoretic force at the high 
salt concentrations used here [22] [22]. 

The flux, J(x , f), is therefore given by a ID convection- 
diffusion equation: 


J( x, t ) = — D 


dc(x, t) 
dx 


c(x, t)n 


dV(x) 

dx 


Dc(x,t) dF(x ) 
kT dx 


(!) 

where D is the diffusion coefficient, x is the centre of mass 
coordinate of the DNA, c(x. t) is the concentration, n is 
the electrophoretic mobility, V{x) is the electric poten¬ 
tial, T is the temperature and F{x) is the entropic bar¬ 
rier. We model the geometry as cylindrical with effective 
length L which includes a small distance in the donor 
compartment where the electric field is strong enough 
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FIG. 3. (a) Modelled ID geometry with length L, concen¬ 

tration Co in the sample reservoir and zero in the opposing 
reservoir, (b) Triangular profile for the entropic barrier F(x) 
with height uo . z; parameterises how far the barrier extends, 
(c) The electric potential (V m . is the applied voltage) and (d) 
combined free energy profile for the nanopore where z is the 
effective DNA charge. 


to trap the DNA (Fig. lb). The entropic barrier is tri¬ 
angular with height uq at a distance r\L/ 2 (Fig. 3). 77 
therefore describes the distance the entropic barrier ex¬ 
tends. We estimated 77=0.25 based on assuming that 
the entropic barrier extends 50 nrn (approximately one 
persistence length of dsDNA) and assuming an effective 
length of L = 200 nm. The concentration is Co at the 
entrance and zero at the exit. 

For single file DNA threading through a nanopore, the¬ 
oretical calculations suggest a weak N~ 0 - 2 dependence 
of barrier height zto on DNA length N [24]. However, 
the nanopore diameter used here allows for non-single 
file threading for which the length dependence of the en¬ 
tropic barrier has not been calculated and we assume a 
barrier height independent of the DNA length. Natu¬ 
rally, a more complicated profile for the nanopore geom¬ 
etry and entropic barrier can be used but we choose this 
one since our goal is to capture the essential physics of 
the transport process and our chosen profile is readily 
solved analytically. 

By applying steady-state boundary conditions, we cal¬ 
culate the translocation frequency as the flux at the 
pore exit from equation ( 1 ). We assume the diffu¬ 
sion coefficient and electrophoretic mobility scale with 
length according to relationships observed in bulk so¬ 
lution. The diffusion coefficient, D , of dsDNA in the 
length regime 0.5 kbp to 10 kbp scales with the chain 
length as D ~ N~ 0 ' 6 [55] [5B]. It is also known that 
electrophoretic mobility is independent of N, /i ~ N° 
122 ■ These dependencies are included in the calculations 
and give the length dependence for the translocation fre¬ 
quency (see Appendix). Our model is appropriate for 
dilute regimes where there are no significant polymer- 
polymer interactions. We performed controls at signifi¬ 
cantly diluted DNA concentrations which showed a simi¬ 
lar dependence of DNA translocation frequency on length 
confirming that we are in this dilute regime (Figures 10- 



FIG. 4. (a) Length dependence of DNA translocation fre¬ 

quency measured at five different voltages (all data from the 
same nanopore) in an electrolyte of 4M LiCl. The trend lines 
show least-squares fits to the data using equation (11) (de¬ 
rived from equation (2) and the model of Fig. 3). (b) Cor¬ 
responding translocation frequency against voltage for four 
selected lengths. Error bars are from the standard errors of 
the Gaussian fits to the ECD histogram. 


12 ). 

Fig. 4a shows our experimental data for the length 
dependence of translocation frequency in 4M LiCl. We 
observe that the translocation frequency increases with 
increasing DNA length with an offset and steepness de¬ 
termined by the applied voltage. The trend lines shown 
are from a least-squares fit to the data in Fig. 4a using 
equation (11). Three parameters; i) the diffusion coeffi¬ 
cient divided by effective pore length, ii) effective DNA 
charge and iii) entropic barrier height were globally fit¬ 
ted to the data set and converge close to estimates from 
literature with an entropic barrier height uq = 4 kT (see 
Appendix). Comparison of the data and fit in Fig. 4a 
shows that the fit deviates for long lengths at low volt¬ 
age and for short lengths at high voltage. However, the 
model accurately captures the trend of increasing translo¬ 
cation frequency with length and the changing gradient 
with voltage. 

Fig. 4b shows the same data as Fig. 4a plotted as 
translocation frequency against voltage for four selected 
lengths. For 0.5 kbp, 2 kbp and 3 kbp the translocation 
frequency is an increasing non-linear function in voltage 
and for 8 kbp it is unclear whether the translocation fre¬ 
quency increases linearly or slightly non-linearly over the 
voltage range. The model fits show a non-linear increase 
for 0.5 kbp, 2 kbp and 3 kbp in agreement with the data 
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and for 8 kbp show a transition from non-linear to lin¬ 
ear though this transition is not clearly observable in the 
data due to the larger error bars for 8 kbp. 

In general the model predicts two regimes of (1) en- 
tropic barrier dominated and (2) drift dominated trans¬ 
port |20j . The entropic barrier regime dominates for 
low voltages and small dsDNA lengths and is charac¬ 
terised by a non-linear increasing translocation frequency 
with voltage and increasing translocation frequency with 
length. The drift regime occurs at high voltage and for 
long dsDNA lengths where equation (1) reduces to give 
a translocation frequency 

J{L) _ nVm 


where V m is the applied voltage and therefore the translo¬ 
cation frequency is linear in voltage and independent of 
dsDNA length. The results presented here in 4M LiCl 
show the characteristics of the entropic barrier regime. 
At higher voltages in Fig. 4a it can be seen that the 
translocation frequency is only slightly increasing with 
length for the longer molecules which suggests that the 
longer DNA lengths at higher voltage are on the thresh¬ 
old of the drift-dominated regime. 

The transition point between barrier-limited and drift- 
dominated regimes is expected to be sensitive to a num¬ 
ber of experimental factors namely pore size and geom¬ 
etry, DNA effective charge (via the type of salt and the 
salt concentration) and temperature. The 4M LiCl elec¬ 
trolyte reduces DNA velocity 10 times compared to 1M 
KC1 (a commonly used electrolyte for nanopores) which 
can be related to a 10 fold decrease in the effective DNA 
charge [25;. We also performed experiments using 1M 
KC1 as the electrolyte where the experimental resolu¬ 
tion is limited to 5 kbp due to the faster translocation 
speed (Fig. 5a). A linear voltage dependence of translo¬ 
cation frequency is observed for 5 kbp, 10 kbp and 20 kbp 
DNA with no significant difference between the translo¬ 
cation frequency for these three lengths indicating a drift- 
dominated behaviour in this length regime and salt con¬ 
centration (Fig. 5b). The lack of a significant barrier 
for translocation can also be observed by extrapolating 
the straight line fit to the data which shows that it in¬ 
tercepts close to the origin. Our results therefore show 
that the threshold voltage for drift-dominated transport 
is significantly decreased in 1M KC1 compared to 4M LiCl 
due to the higher effective charge. These results at 1M 
KC1 are supported by experiments with similar diameter 
nanopores which showed a linear dependence for A-DNA 
(48.5 kbp) translocation frequency on voltage at 1M KC1 
[29] . Also results for a pressure-voltage trap system, with 
comparable conditions of 1.6M KC1 and a 10 nm diameter 
pore, show that there are expected to be very few unsuc¬ 
cessful translocation attempts for DNA lengths of a few 
kilobase-pair under voltage-only experiments and hence 
that the electrophoretic drift dominates at this regime 

m- 
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FIG. 5. (a) Scatter plot of 1168 translocations from a sam¬ 

ple containing 5 kbp, 10 kbp and 20 kbp with each length 
at a concentration of 0.12 nM and using 1M KC1 electrolyte 
with +600 mV applied voltage. The dashed grey lines are 
fits according to A I Mean = — ECI ^' lak and the solid line 
shows the 3 fC threshold, (b) Corresponding translocation 
frequency against voltage for 5 kbp, 10 kbp and 20 kbp - the 
translocation frequency values are significantly higher than 
that observed in 4M LiCl (Fig. 4b) due to the higher effec¬ 
tive DNA charge in 1M KC1. The line shows a linear fit to 
all data points yielding a gradient of 0.011 nM" 1 s~ 1 mV~ 1 . 
Error bars are from the standard errors of the Gaussian fits 
to the ECD histogram. 


V. CONCLUSIONS 

To conclude, we have investigated the translocation fre¬ 
quency of dsDNA through quartz nanopores with mean 
diameter 15 nm. The low noise of these nanopores, com¬ 
bined with the narrow distribution of dwell times for 
individual lengths in such wide pores, allows us to ac¬ 
curately measure the translocation frequency using cal¬ 
ibrated DNA samples. At an electrolyte concentration 
of 4M LiCl we have observed an increasing translocation 
frequency as a function of DNA length in the range 0.5 
kbp to 10 kbp which fits the behaviour of entropic barrier- 
limited transport in a ID convection-diffusion equation 
and allows us to extract an approximate value of 4 kT for 
the entropic barrier in this model. At a lower salt con¬ 
centration of 1M KC1 where the effective DNA charge 
is higher, we measure the characteristic transport prop¬ 
erties of a drift-dominated regime. Our model accounts 
for the behaviour in these two different regimes and pro¬ 
vides a basis for understanding and predicting polymer 
transport through nanopores across a wide range of con- 
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APPENDIX: ANALYTICAL MODEL 
DERIVATION 


Equation (1) gives the ID convection-diffusion equation: 


T , . ^dcix^t) . . dVix) 

J(x, t) = —D — - c(x, t)iu- 


dx 


dx 


the mobility p. is given by 


zD „„ 

"= kf =iD ' 


Dc{x, t) dF(x) 
kT dx 


( 3 ) 


( 4 ) 


where z is the effective charge and z = -Sp . Also writing 
F(x) = ^r- gives an expression for the flux of 


J(x, t ) = — D 


dc(x,t ) . dV(x) . .dF(x) 

v ; +c{x,t)z—^ +c{x,t) 


dx 


dx 


( 5 ) 

To calculate the translocation frequency we find the 
steady state solution of this equation with the bound¬ 
ary conditions c(x = 0) = eg, c(x = L) = 0. The en- 
tropic barrier F(x) is assumed to be triangular as shown 
in Figure 3 so that: 


j 2u 0 f , 0 < x < 4 

F(x) = l 2uo(l - f), | < x < l (6) 

[ 0, l < x < L 

where l = rjL and uq = The electrical potential in 
the pore is given by 

V(x) = V m 0 <x<L (7) 

where V m is the applied potential. Therefore the total 
free energy divided by kT, W(x), is given by 

_ f 2 uof — \zV m \ x, 0<x<4 
W(x) = l 255(1 - f) - \zV m \ f, ±<x<l (8) 

( -\zV m \f;, l < X < L 

In steady-state the flux is independent of time and x and 
is calculated by integrating equation (5) which gives [2D] 


J{L) = -D 


c(L)e 


W(L) _ 


c( 0 )e 


W( 0) 


f 0 e w ( x )dx 


( 9 ) 



FIG. 6. Diffusion coefficient of DNA against N ~°' 6 where N 
is in basepairs (using data from |25|). A least-squares linear 
fit to the data is shown. 


and since c(x = 0) = cq, c(x = L) = 0, 


J(L) = D 


co 


f 0 e w ( x )dx 


( 10 ) 


calculating this integral with the assumed free energy 
profile W(x) gives 




L ' \(^ !1 - \zVm\) 




e v 77 


- 1 ] 


(** + \ZV m \) 


^ + \JVm\)r, _ e -{=^M-|SV m |)§i 


\w m 


[ e -|2V m | _ e -\zV m \m 


-1 


( 11 ) 


When the electrophoretic force dominates over the bar¬ 
rier 77 \ zV m \ Mo and equation ( 11 ) reduces to 


J(L) = 


Dc 0 zV m CQ^Vm 


( 12 ) 


which is the same as equation ( 2 ) and is the drift- 
dominated regime found in 1M KC1 in Figure 5. When 
the drift does not dominate, we must use the full equation 
(11). The length dependence enters from the dependence 
on D and z, and we consider these in turn. The diffusion 
coefficient scales according to [25) 25 , 

D = D 0 N~ 0 - 6 . (13) 


Sorlie et al. used dynamic light scattering to determine 
the diffusion coefficient of isolated DNA molecules in 
100 mM NaCl, 10 mM Tris-HCl (pH= 8 ), 1 mM EDTA 
at 20°C. Their data yields a value of D 0 = 5.9 x 
10 -10 m 2 s~ 1 (Figure 6 ). 

The electrophoretic mobility of double-stranded DNA 
> 100 bp (and other polyelectrolytes) is experimentally 
measured to be independent of length [27] I30H32] 

d ~ N°, 


(14) 
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and therefore 

z = z 0 N 0 - 6 . (15) 

We estimate Sq by the following argument. DNA has a 
bare charge from its phosphate groups of 2e“/bp. Exper¬ 
iments on tethered DNA in a nanopore at 1M KC1 showed 
a 75% reduction in charge comparable to the charge re¬ 
duction expected from Manning condensation [33]. Fur¬ 
thermore, as reported by Kowalczyk et al. [28], we ob¬ 
serve a ten-fold reduction in DNA velocity for experi¬ 
ments in 4M LiCl compared to 1M KC1 for our nanopores. 
Therefore we estimate an additional factor of 10 reduc¬ 
tion in effective charge giving a total of a 40 fold re¬ 
duction from the bare charge. Experiments are done at 
20°C so kT = 4.0 x 10~ 21 J and the total charge on 2e“ 
is 3.2 x 10 -19 C. Therefore we obtain So = 2 CJ~ l . This 
value of So together with that estimated for Do yields a 
mobility of ^ = SqDq = 1.2 x 10 -9 m 2 s~ 1 V~ 1 which 
is comparable to a value of /j, = 1.8 x 10~ 8 ?n 2 s~ 1 V~ 1 
determined for 10 kbp DNA in 1M KC1 01] given the 
anticipated ten fold reduction due to 4M LiCl. 

The other parameter in our model is the barrier height. 
Previously the barrier height was estimated for single 
chain threading to be on the order 10-20 kT |24]. The 
pores here allow for folded configurations as the DNA 
passes through and to reflect this we use initial estimates 
of 3-10 kT for the barrier height. Furthermore, as de¬ 
scribed in Section IV, we assume that the barrier height 
is independent of the DNA length. 

We fit equation (11) to the data in 4M LiCl using three 
parameters: a — D 0 /L (where L is the effective length 
taken as L = 200 nm ), So , and uq using the initial values 
a = 0.003 ms- 1 , So = 2 CJ -1 and uq = 3 — 10 in steps 
of 1. A least-squares global fit, using the Levenberg- 
Marquadt algorithm (OriginPro 8.5), to Figure 4a yields 
a = 0.061 ms -1 , So = 0.58 CJ -1 and Uq = 4.3 which 
converged for all values of tested values of Uq. The pa¬ 
rameter a changes most significantly relative to the ini¬ 
tial guesses. We note that this change can be expected 
given the uncertainty in the effective pore length as well 
as assumptions about other parameters such as the bar¬ 
rier width, shape and length dependence which were es¬ 
timated as described in the text. Also a larger diffusion 
coefficient (and therefore a) might be expected compared 
to the initial guess since we calculated D 0 based on ex¬ 
periments at 100 mM NaCl by Sorlie et al. whereas the 
nanopore experiments presented here use 4 M LiCl elec¬ 
trolyte and high salt concentrations can cause an increase 
in the polymer diffusion coefficient [35]. We would like to 
note that our model requires a length dependent diffusion 
constant to describe the data for the DNA translocation 
frequency at 4M LiCl since /i and F(x) are length inde¬ 
pendent. If D were also length independent, equation (1) 
would then have no N dependence which is contrary to 
our observations. 

APPENDIX: FURTHER EXPERIMENTAL 
CHARACTERISATION 



- 0.6 - 0.4 - 0.2 0.0 

Current change (nA) 


FIG. 7. (a) Typical translocation events of 10 kbp DNA show¬ 
ing steps due to the folding configurations of the DNA as it 
translocates the nanopore, (b) All points histogram of 260 
events showing a baseline peak at 0 nA and successive peaks 
due to DNA folding states. 


Nanopore geometry 


Nanopores in this study were fabricated with the same 
protocol as given in Bell et al. m where the pore diam¬ 
eters were estimated using scanning electron microscopy 
as 15 ± 3 nm (mean ± sd). We can also obtain an indica¬ 
tion for the pore size by analysing the number of folded 
DNA translocations. Figure 7 shows translocation events 
and an all points histogram from 10 kbp translocations 
at +600 mV (selected by ECD from the data shown in 
Figure 2). There is a peak at 0 pA due to the the base¬ 
line level and the one DNA double-strand level occurs at 
-122 pA. A second peak is also visible at -217 pA due 
to configurations where the DNA has a fold so that two 
double-strands are present within the pore. A third peak 
is also visible confirming that the pore diameter allows 
at least three double-strands and is therefore consistent 
with a diameter of ~15 nm. 











600 



0 100 200 300 

Event charge deficit (fC) 


FIG. 8. Example least-squares fit of multiple Gaussian func¬ 
tions to the 10 peaks from the 10 DNA lengths present in the 
sample given in Table I. The data is the same as that present 
in Figure 2 but with a linear rather than logarthmic binning. 


Determination of number of translocations from 
ECD histogram and time between successive events 

The number of DNA strands that translocated in an ex¬ 
periment was determined by fitting multiple Gaussian 
functions to a histogram of ECD values. The fitting was 
done using the Levenberg-Marquadt algorithm to deter¬ 
mine a least squares fit to the distribution (OriginPro 
8.5). The number of peaks and approximate peak posi¬ 
tions were given as initial parameters. Figure 8 shows an 
example fit to an ECD histogram data set. The ECD 
range was limited between the threshold of 3 fC and 
/i+4cr for the final peak. Figure 9 shows the time between 
successive events for the 5711 events presented in Figure 
8 . We observe a clear exponential decay as expected for 
a Poisson process with independent DNA translocations 
such that the probability distribution function for a ar¬ 
rival in time t after a translocation is 

P(t) oc e~ xt . (16) 


Dependence of DNA translocation frequency on 
concentration 

We tested the translocation frequency dependence at dif¬ 
ferent concentrations relative to those given in Table I. 
Figure 10 shows the translocation frequency at the con¬ 
centration given in Table I (RC=1), and at three times 
dilution (RC=0.33). These concentrations were tested 
using the same nanopore and at +600 mV applied volt¬ 
age. We observe no significant difference after dilution of 
the DNA sample. 

We also performed experiments on a custom-made 
DNA ladder to check for potential sample variability and 
to further extend the concentration range tested. The 



Time between successive events (s) 


FIG. 9. Histogram of time since last event for all 5711 events 
shown in Figure 2. The solid line is a fit according to equation 
16. 



FIG. 10. Dependence of DNA translocation frequency on 
DNA concentration. Translocation frequency against DNA 
length measured for an individual nanopore at the concentra¬ 
tion given in Table I (RC=1) and at a 3 times dilution relative 
to Table I (RC=0.33). 


ladder was made by mixing 8 chromatography-purified, 
dsDNA lengths (NoLimits, ThermoScientific, UK). Each 
DNA length concentration was measured from the ab¬ 
sorbance at 260 nm using a Nanodrop spectrophotome¬ 
ter. The DNA lengths were then mixed together and 
diluted in lx TE, 4 M LiCl so that the final concentra¬ 
tion of each DNA length in the sample reservoir was as 
given in Table II. Figure 11 shows an ECD histogram 
and scatter plot of 1066 events taken in 21 minutes of 
recording at +600 mV with this DNA sample. Figure 12 
shows the translocation frequency of DNA correspond¬ 
ing to the data set in Figure 11 and an additional data 
set from another nanopore. We measure the same trend 
as Figure 4 of an increasing translocation frequency with 
DNA length. The concentration was significantly lower 
than that used in Table I and resulted in ~1 translo¬ 
cation per second making it highly unlikely that there 
were effects due to interactions of DNA strands at the 
nanopore mouth. As an experimental note, we did not 
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DNA 

length (kbp) 

1 

2 

3 

5 

7 

10 

15 

20 

Cone (nM) 

0.22 

0.21 

0.21 

0.19 

0.20 

0.18 

0.18 

0.18 


TABLE II. Final reservoir concentration for DNA lengths in 
custom DNA ladder used for Figures 11 and 12. 


(a) 



Event duration (ms) 

( b ) 20kbp 



Event charge deficit (fC) 


FIG. 11. Experiment on custom ladder given in Table II at 
4M LiCl. (a) Scatter plot of all 1066 translocations. Dashed 
lines indicate constant ECD with values determined from the 
centre of Gaussian peaks fitted to each ECD histogram, (b) 
ECD histogram of 1066 translocation events of the DNA lad¬ 
der given in Table II. Eight peaks can be distinguished corre¬ 
sponding to the eight DNA lengths present in the sample. 

use lengths longer than 20 kbp in this study (such as 
48.5 kbp A-DNA) due to the difficulty in preparing accu¬ 
rate concentrations of such lengths. This stems from two 
sources 1) an increased amount of fragments is observed 
with such long lengths and 2) the slow rate of mixing 
due to the small diffusion coefficient which means that 
inhomogeneous concentrations can easily arise. 

Repeats of voltage dependence using different 
nanopores 

Figures 13 and 14 shows repeats, using different 
nanopores, of the voltage dependence behaviour shown in 



FIG. 12. Translocation frequency as a function of DNA length 
for the DNA sample given in Table II at +600 mV. The two 
colours represent data from two nanopores with fits yielding 
a = 0.10 ms _1 , So = 0.38 CJ _1 and uo = 4.1 (blue data) and 
a = 0.12 ms~ 1 ,So = 0.31 CJ _1 and Uq = 4.6 (brown data). 
Error bars represent the standard errors of the fits. 

Figure 4 (using the DNA sample given in Table I). Each 
nanopore can be expected to give slightly different val¬ 
ues for translocation frequency due to small differences in 
geometry. However we consistently observe, in 4M LiCl, 
the trend of increasing translocation frequency with in¬ 
creasing DNA length. We also consistently observe that 
the data for 4 kbp is systematically slightly above the 
general trend which could be due to a slight error in the 
DNA concentrations or due to the fact that the 4 kbp 
ECD peak is next to the large 3 kbp peak originating 
from the high 3 kbp concentration. A slight deviation 
from a Gaussian peak could therefore affect the calcu¬ 
lated 4 kbp frequency. 

Translocation statistics lie oustide analysis 
thresholds 

The thresholds of detection of 50 pA deviation from base¬ 
line current, minimum 50 /is event duration and min¬ 
imum 3 fC event charge deficit were used in all experi¬ 
ments. We carefully checked for each experiment that the 
translocations lied outside these thresholds and therefore 
that there are no missed translocations which would cre¬ 
ate systematic error in the calculated translocation fre¬ 
quency. For example, Figure 15 shows scatter plots of 
mean event current and peak event current against event 
duration for the particular nanopore used in Figure 2 
and 4. The 50 pA and 50 /is thresholds are marked for 
the peak current scatter plots and the 3 fC threshold is 
marked for the mean event current scatter plots. The 
main clusters of points for different lengths lie outside 
these limits therefore indicating that the DNA lengths 
down to 500 bp are accurately measured. 
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FIG. 13. Repeat of Figure 4 with a different nanopore show¬ 
ing the consistent behaviour of increasing translocation fre¬ 
quency with increasing DNA length. The global least-squares 
fit yields a = 0.073 ms -1 , zo = 0.68 CJ _1 , uq = 4.5. 


FIG. 14. Additional repeat of Figure 4 with another nanopore 
with fits yielding a = 0.055 ms -1 , zq = 0.56 CJ _1 , uo = 4.3. 
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FIG. 15. Left column: Mean event current against event 
duration for all translocations. The green line shows the 3 fC 
threshold set in the analysis. Dashed lines indicate constant 
ECD with values determined from the centre of Gaussian 
peaks fitted to each ECD histogram. N=637 (+300 mV), 
N=3174 (+500 mV) and N=11661 (+700 mV). +300 mV and 
+500 mV were recorded for approximately 20 minutes and 
+700 mV was recorded for approximately 30 mins. Right 
column: Corresponding plots of peak current value against 
event duration. The threshold limits of 50 fis and 50 pA are 
also marked as green lines. Events cluster outside the thresh¬ 
old limits indicating that nearly no translocations are missed. 
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